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ABSTRACT 

Highly magnetized neutron stars (NSs) are characterized by a bewildering range of astrophysical manifesta- 
tions. Here, building on our simulations of the evolution of magnetic stresses in the NS crust and its ensuing 
fractures (Perna & Pons 2011), we explore in detail, for the middle-age and old NSs, the dependence of star- 
quake frequency and energetics on the relative strength of the poloidal (Bp) and toroidal (Btoi) components. We 
find that, for Bp > lO''' G, since a strong crustal toroidal field Btoi ^ Bp is quickly formed on a Hall timescale, 
the initial toroidal field needs to be Btoi ^ Bp to have a clear influence on the outbursting behaviour. For initial 
fields Bp < 10''* G, it is very unlikely that a middle-age (f 10^ years) NS shows any bursting activity. This 
study allows us to solve the apparent puzzle of how NSs with similar dipolar magnetic fields can behave in a 
remarkably different way: an outbursting 'magnetar' with a high X-ray luminosity, or a quiet, low-luminosity, 
'high-B" radio pulsar. As an example, we consider the specific cases of the magnetar IE 2259H-586 and the 
radiopulsarPSR J18 14-1744, which at present have a similar dipolar field ^ 6 x 10'^ G. We determine for each 
object an initial magnetic field configuration that reproduces the observed timing parameters at their current 
age. The same two configurations also account for the differences in quiescent X-ray luminosity and for the 
'magnetar/outbursting' behaviour of IE 2259H-586 but not of PSR J1814-1744. We further use the theoret- 
ically predicted surface temperature distribution to compute the light-curve for these objects. In the case of 
IE 2259H-586, for which data are available, our predicted temperature distribution gives rise to a pulse profile 
whose double-peaked nature and modulation level is consistent with the observations. 

Subject headings: stars: magnetars — stars: neutron — X-rays: bursts 



1. INTRODUCTION 

Neutron stars, end points of the evolution of massive stars, 
are characterized by a bewildering variety of astrophysical 
manifestations. Among these, particularly intriguing is a class 
of slow rotators (P ^ 2 - 1 1 sec) with an especially high qui- 
escent X-ray luminosity (L^ ^ lO^-'- 10^^ erg), almost always 
exceeding their rotational energy loss E, and displaying oc- 
casional energetic outbursts and giant flares. The latter can 
release energies as large as 10"*^ erg. 

The source of energy responsible for powering the out- 
bursts and flares, and enhancing the quiescent emission, is be- 
lieved to be of magnetic origin (Thompson & Duncan 1995, 
2001). The magnetic field evolves in the NS interior un- 
der the effects of Ohmic dissipation. Hall drift and Ambipo- 
lar diffusion, although this latter has recently been suggested 
to play a lesser role in the presence of superfluid matter 
(iGlampedakis. Jon es & Samuelsson 2011). In this case, the 
evolution of the magnetic field causing the magnetar activity 
would likely take place in the stellar crust. However, if the ob- 
servational evidence from the cooling of the supernova re m- 
nant in Cassiopeia A (Pag e et al.l201 ItlShternin et al.l201 ih is 
confirmed, it would imply that the object is undergoing a core 
superfluid transition at 5 x 10^ K. In the magnetar context, 
higher core temperatures are expected and the transition to 
core superflu idity can be significantly delayed, a s originally 
suggested by Arras. Gumming & ThompsonI (1200 4). There- 
fore, the role of ambipolar diffusion in the core remains an 
open question. In this work we focus on the role of Ohmic 
diffusion and Hall drift in the crust. 
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Occasionally the magnetic stress exceeds locally the tensile 
strength of the crust, causing it to fracture; the ensuing sud- 
den release of elastic and magnetic energy powers the out- 
bursts. This picture, in which the magnetic field strength 
'regulates' the properties of an object, i.e. its appearance 
as an outbursting magnetar, or a 'normal', quiet radio pul- 
sar, left many unanswered questions and unsolved puzzles. 
Why some objects are more active than others? What deter- 
mines the frequency and energetics of events? Why some ob- 
jects have been observed to emit very energetic 7-ray flares 
(the so-called "Soft Gamma-ray Repeater", (SGRs)), while 
others only display less energetic X-ray bursts (the so-called 
"Anomalous X-ray Pulsars" (AXPs))? Is there a real differ- 
ence between these two classes ? And, if the magnetic field 
is the main driving force, how is it possible for a "low B"- 
field pulsar to undergo an outburst (Rea et al. 2010)? Alto- 
gether, how can two objects, characterized by an almost iden- 
tical magnetic field (as inferred through measurements of the 
pulsar period P and its derivative P) manifest themselves as 
apparently different objects with respect to those properties 
which should be driven by the B-field strength (outbursting 
behaviour, quiescent X-ray luminosity)? This apparent di- 
chotomy is exemplified by the magnetar IE 2259-1-586, with 
an infeiTed dipolai" field of 5.9 x 10'^ G, and the 'high-B' 
field radio pulsar PSR J1814-1744, with an inferred Bdip = 
5.5 X lO'^* G. The former is classified as an AXP; it was ob- 
served to have an active period of outbursts in 2002 (Woods 
et al. 2004), and has an enhanced X-ray luminosity in qui- 
escence, Li ^ lO^"* erg, which exceeds its rotational energy 
loss. On the other hand, PSR J1814-1744 appears to be a 
"normal", quiet radio pulsar (Camilo et al. 2000), with no 
discernible X-ray emission in the X-ray band. 

In order to understand from first principles what is the phys- 
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Fig. 1 . — Luminosity at two different ages as a function of tlie initial magnetic field. Left panel: dependence of the luminosity on the initial poloidal field (at 
pole) for models with Btoi = 0. Right panel: same for models with Bp = 3 X lO''' G and different strength of the initial toroidal field. 



ical origin of the variety of observed phenomenology, in Perna 
& Pons (2011; PP2011 in the following), we performed the 
first long-term 2D simulations that followed the evolution of 
magnetic stresses in the NS crust; these, combined with cal- 
culations of the breaking stress of the NS, allowed us to theo- 
retically predict frequency, energetics, and location on the NS 
surface of the outbursts (associated with crust fractures). Our 
results allowed us to establish that there is no fundamental dif- 
ference among apparently different manifestations of the ob- 
jects (such as 'AXPs' and 'SGRs'). We also found that both 
the dipolar component and the 'hidden' toroidal component 
play an equally important role in determining the frequency 
and energetics of the starquakes. Furthermore, for a given ini- 
tial B-field configuration and strength, younger objects yield 
more frequent and more energetic outbursts. 

In this follow up paper, we continue our work in several 
respects. Firstly, we extend our study of the outburst fre- 
quency and energetics to focus on the influence of the rela- 
tive strengths of the poloidal and toroidal field components. 
This study allows us to uncover some fundamental clues as 
to what creates the varied phenomenology of the highly mag- 
netized neutron star family. Second, for the middle-aged and 
older objects, we consider the mutual influence of the mag- 
netic and thermal evolution and predict, for a combination of 
Bp and Btm initial strengths, the X-ray luminosity, tempera- 
ture, timing properties and magnetic field at current age. Last, 
as an example, we consider two specific objects, IE 2259-1-586 
and PSR J1814-1744, which, as described above, have almost 
identical spin parameters and hence inferred dipolar B field, 
but very different manifestations. For each of them, we show 
how we can find an initial B-field configuration that repro- 
duces their observed timing parameters at their current timing 
age, as well as account for the differences in quiescent X- 
ray luminosity, and for the 'magnetar/outbursting' behaviour 
of IE 2259H-586 but not PSR J1814-1744. We also compute 



the surface temperature distributions associated with the mag- 
netic configurations that reproduce those properties and use 
them to compute the expected X-ray light-curves in quies- 
cence. For the case of IE 2259-1-586, for which data are avail- 
able, we show that the qualitative features that we predict are 
a good match to the data. 

The paper is organized as follows. The magneto-thermal 
evolution is described in §2, together with predictions for the 
intensity of the X-ray luminosity with age as a function of 
Bp and Btor- In §3, we study the dependence of the starquake 
frequency and energetics on the relative strengths of Bp and 
Btoi . The specific examples of magneto-thermal evolution for 
IE 2259H-586 andPSR J1814-1744 are presented in §4. We 
summarize and discuss our work in §5. 

2. MAGNETO-THERMAL EVOLUTION OF HIGHLY MAGNETIZED 
NEUTRON STARS 

The first attempt to study the long term cooling of mag- 
netized NSs, including all realistic microphy sics and taking 
into account magnetic field decay was made bv lAguilera et alj 
(2008). They assumed a simple, analytical law for the time 
variation of Bp which incorporates the main features of the 
Ohmic and Hall terms in the induction equation. Numer- 
ical simulations of the magnetic field evolution in a real- 
istic neutron star crust, including the nonlinear effec ts of 
the Hall term were presented in iPons & Gepperj (12007 ^. but 
in this case assuming an isothermal crust. The fully cou- 
pled magneto-thermal evo lution of a NS was later studied by 
iPons. Miralles & Gepperll (2009) . However, owing to numer- 
ical difficulties in treating the Hall term, their models could 
include only Ohmic diffusion. A fully coupled evolution- 
ary simulation has not been possible up to now, but for ini- 
tial values of Bp < 10'^ G, the effect of the Hafl drift is 
expected to introduce at most quantitative changes (a some- 
what faster dissipation) with respect to the purely Ohmic pic- 
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Fig. 2. — Average number of outbursts per 1000 years as a function of 
the initial toroidal field strength, for a lO' years old neutron star. We show 
results for thi'ee different initial poloidal fields: Bp(t = 0) = lO" G (soHd 
lines), Bp(t = 0) = 5 X 10"' G (dashes), and = 0) = 3 X lO'-* G (dash- 
dotted lines). 



ture described in lPons. Miralles & Gepper3(l2009h . which can 
be considered a good approximation to reality. With this 
code, we performed a number of simulations of our base- 
line model (a IAMq NS with a radius of 11.6 km, and mo- 
ment of inertia / = 1.4 x 10^^ g cm-) varying the initial val- 
ues of the dipolar and in t ernal t oroid al fields. We refer to 
iPons. Miralles & G epperll (|2009|) and lAguilera et all (l2008h 
for all technical details about the cod e and the microp hysical 
input. The minimal cooling scenario ('P age et ani2004h . with- 
out exotic phases nor fast neutrino cooling processes, was as- 
sumed. We also include enhanced neutrino emission from the 
breaking and formation of neutron Cooper pairs in the core, 
consistently with the observational evidence from the cooling 
of the supernova rem nant in Cassiopeia A dPage et al.ll201 ll : 
IShternin etal]l2011h . 

As noted in previous studies, the average luminosity of a NS 
born with Bp < 5 x 10'^ G. is hardly affected, compared to the 
non-magnetized case. On the other hand, those NSs born as 
magnetars are subject to significant heating by the dissipation 
of currents in the crust. At fixed poloidal field, the luminosity 
of models with strong internal toroidal components is system- 
atically higher than that of models without toroidal fields, due 
to the additional energy reservoir stored in the toroidal field, 
which is gradually released as the field dissipates. In Fig. [T] 
we show the dependence of the luminosity at a fixed age as 
a function of the initial magnetic field. We can see that, for 
low initial dipolar fields, the luminosity barely changes with 
increasing dipolar field but, once we enter in the magnetar 
range, Bp > 10'^ G, a further increase in Bp results in a lumi- 
nosity increase of more than one order of magnitude. At fixed 
poloidal field, the existence of an internal toroidal field has no 
evident effect if Btor ^ Bp, but the luminosity is again largely 
increased when the internal toroidal field exceeds Bp. This 
can be understood with the following back-of-the-envelope 
estimate. The rate at which magnetic energy e„, is dissipated 
can be approximated by 
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where J = xB is the current density, V is the volume 
where currents are dissipated (the crust volume) and a is the 
electrical conductivity. Approximating V x B « B/L, where 
L is a typical scale-length, and for typical neutron star condi- 
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Fig. 3. — Outburst energy distribution of objects at age lO' yr bom with 
Bp = 5 X lO'"* G and four different initial toroidal fields: Btor = (blue dots), 
10'^ G (red dashes), lO"^ G (green dash-dotted line), and 3 X lO"^ G (black 
solid line). The total number of outbursts corresponds to a period of 1000 yr. 



tions we have, for L = 1 km, that 

de, 



1033^erg/s. 

dt (724 



(2) 



Here B15 is the magnetic field in units of 10^^ G and (T24 is 
the electrical conductivity in units of lO'** s"'. In order for 
the internal field to have a significant effect on the luminosity, 
the dissipation rate of magnetic energy has to be larger than 
the luminosity in the non-magnetic case. This explains the 
right panel of Fig. [T] where the influence of the toroidal field 
only becomes visible after a critical value of 10'^ G and then 
increases rapidly with the field strength. 

We must remind that the quoted values of Bp and B^f corre- 
spond to the initial values. In the models with stronger fields, 
as a result of the average higher temperatures, the crustal elec- 
trical resistivity is also enhanced and magnetic diffusion is 
more efficient during the first 10^-10'' years of a NS's life. 
As a consequence, all NSs born with fields > 5 x 10''* G 
are subject to faster dissipation. Typically, the dipolar field 
decreases by a factor of 2-3 in the first 10^ years of a NS 
life, and it asymptotically reaches a similar field strength 
(« 2-3 X lO''' G) at very late times, irrespectively of the ini- 
tial strength. The toroidal field at 10'' years is reduced (on 
average) by one order of magnitude with respect to the initial 
value. Note also that, even if the initial model has no toroidal 
field, for sufficiently large (poloidal fields > 5 x 10'"* G) the 
Hall drift will create a toroidal field of similar strength. 

3. OUTBURST FREQUENCIES AS A FUNCTION OF THE STRENGTH 
OF THE DIPOLAR AND TOROIDAL B-FIELD COMPONENTS 



In P P201 1, we used the numerical code bv lPons & Gepperl 
(|2007|) to follow in axial symmetry the evolution of the mag- 
netic field in the crust of the neutron star. The magnetic stress 
Mij{r,9,t) = Bi(r,9,t)Bj(r,9,t)/4Tr is computed at each time 
as a function of the latitude 6 on the surface of the star, and 
the radial coordinate r in the NS crust. As the magnetic field 
evolves, the crust moves through a series of equilibrium states, 
during which elastic stresses balance variations in magnetic 
stresses. If, during the evolution, any component of the mag- 
netic stress departs from its equilibrium value by an amount 
comparable (or exceeding) the breaking stress of the crust, the 
crust will break, releasing the accumulated elastic energy (for 
technical details of the calculations, see PP2011). The simu- 
lations allowed PP201 1 to predict the frequency and energet- 
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ics of the starquakes, for a certain initial Z?-field configuration 
characterized by a poloidal and a toroidal component. The 
case of Bp = 8 X lO'^f G and Btoi- = 2 x 10'^ G was considered in 
detail in that work. It was shown how both the frequency and 
the energetic of the outbursts decrease with the NS age. Some 
initial exploration of the role of different initial field strengths 
further showed that the toroidal component does play a very 
important role in the outburst properties. 

The former study was limited by the numerical restrictions 
of the code that do not allow to perform a systematic study 
of the starquake frequency and energetics of objects with 
> 10^ yr with the non-linear terms (too long runs). In or- 
der to overcome the numerical difficulties we have adopted 
a hybrid approach. Since, for these objects, the Hall time 
is much shorter than their age, we expect that the initial 
configuration relaxes to a nearly Hall-neutral geometry on 
a few Hall timescales (reaii = A-nrieeL^ /(cB) ^ 6000B7{ yr), 
after w hich the evolu tion becomes dominated by Ohmic dif- 
fusion (iPons & Gepp ert 2007). Thus, as a first approxima- 
tion, we follow the NS evolution up to 10^ yr with the code 
of Pons, Miralles & Geppert (2009), and then use the temper- 
ature and magnetic field profiles as initial input for the short- 
term, fully non-linear e volution including the Ha ll term, using 
the numerical code by iPons & Geppert! (l2007h . The caveat 
of this approach is that it does not account for possible fast. 
Hall-induced, initial changes of the magnetic field geometry 
and strength that may occur for magnetar conditions. This can 
result in a dissipa tion of the initial field in , at most, a factor of 
2, as estimated in lPons & Gepperj(l2007l) . It is also important 
to note that, whatever the initial geometry of the toroidal field 
was, it is rapidly (Hall timescale) driven to a similar shape 
that corresponds to a sort of quasi-steady state. Therefore, 
our results for middle aged NSs (> 10^ yr) are expected to be 
rather insensitive to the particular form of the toroidal field, 
being more relevant the role of global quantities (total energy 
or helicity) that are explored by varying the normalization of 
the toroidal component. 

Our main results concerning the outburst frequency and en- 
ergetics are summarized in Fig. |2] and Fig. 
bfreffig3. In Fig. |2] we show the dependence of the average 
number of expected outbursts in a period of 1000 yr on the 
strength of the initial magnetic field. We have taken the pa- 
rameter e = 0.9 (see PP201 1). We can see that the average out- 
burst rate is sensitive to both Bp and Btor in a similar way as the 
luminosity. Below a critical poloidal field Bp < 3 x 10'"* G the 
burst rate is vanishingly small (< 10"^ yr~'), except for very 
large toroidal fields Btor > 10'* G, but even in that extreme 
case the outburst frequency is smaller than one per century. 
For initial Bp > 3 x 10'"*, the outburst rate grows exponen- 
tially with the value of Bp. At a fixed Bp, the dependence on 
the toroidal field is very weak until Btor ^ Bp, but the out- 
burst rate again increases rapidly with increasing values of 
Btor > Bp. It should be noted that these numbers must be taken 
with caution, simply as an order of magnitude estimate. The 
large uncertainties about the internal geometry of the mag- 
netic field and details of the crustal fracture mechanisms do 
not allow a precise calculation. Nevertheless, these estimates 
may serve as an indication of the average outburst frequency 
for objects with similar conditions. The energy distribution 
of the events during a period spanning 1000 years is shown 
in Fig. [3] for the model with Bp = 5 x lO''* G and four dif- 
ferent toroidal fields. Except for the model with the largest 
Btoi, there is no significant variation in the outburst average 



event energy, which is predominantly in the range 10"*' - 10"^^ 
erg. Only for the extreme case of B^r = 3 x 10'* G, in addi- 
tion to the large increase in the number of events, the average 
outburst energy is also substantially larger, up to 10'*-' erg, re- 
flecting the larger energy reservoir available in the toroidal 
field. 

4. IE 2259+586 AND PSR J1814-1744: SIMILAR Bp BUT DIFFERENT 
ASTROPHYSICAL MANIFESTATION 

The comprehensive study performed in the previous two 
sections allows us to obtain a physical explanation to an ap- 
parent puzzle in neutron star theory. How can two objects 
have a very similar dipolar field but display a markedly differ- 
ent behaviour, if the magnetic field provides the main driving 
force? 

Our magnetothermal simulations provide an answer to this 
question: for a similar inferred value of the dipolar compo- 
nent of the magnetic field, the strength of the internal toroidal 
field can have a dramatic influence on the behaviour of the 
neutron star, and in particular on its manifestation as an X- 
ray bright, outbursting 'magnetar', or an X-ray dim, quiet 
'high-B' radio pulsar In the previous section we have shown 
that, if the initial dipolar component is low (< x lO''* G), the 
toroidal component, even if strong, is not very effective in 
producing frequent outbursts. On the other hand, for dipo- 
lar field strengths on the order of a few x 10'^ G, a toroidal 
component of the same strength is produced during the evo- 
lution, even if it was not there initially. Thus, the impor- 
tance of an existing initial Btor is only evident when it ex- 
ceeds the poloidal component. As our results have shown, 
over ages of about 10^ yr, a dipolar B-field of a few 10'** G 
roughly halves. It is especially interesting to note that it is 
around dipolar magnetic fields of 6-7 x 10'^ G (and ages of 
about 10 yr) that objects classified as 'magnetars' or 'high- 
B' radio pulsars are found to coexist. The most notable ex- 
amples, as discussed in §1, are IE 2259H-586, with an in- 
ferred Bp = 5.9 X 10'^ G and magnetar characteristics, and 
PSR J1814-1744, with Bp = 5.5 x 10'^ G and properties of a 
"normal" radio pulsar We do not include other interesting ob- 
jects in this study, such as the high magnetic field radio pulsars 
PSR J1119-6127 or PSR J1846-0258, because they ai-e very 
young (w 1 kyr) and a fully coupled treatment of the Hall term 
with the thermal evolution would be required. We have how- 
ever also examined the high-B radio pulsar PSR J1718-3718, 
and found that the properties of this object can be well ex- 
plained by the same model used for PSR J1814-1744, but in 
an earlier stage. In the remaining of this section we specif- 
ically focus on the cases of IE 2259H-586and PSR J1814- 
1744and we leave a more extensive study of many different 
objects for the future. 

We evolved in time our baseline NS varying only the initial 
B-field configurations (characterized by Bp and B^r) trying to 
reproduce the timing, the magnetic, and the thermal properties 
of these two objects at their measured timing age. We assume 
an orthogonal rotator (angle between the magnetic axis and 
rotation axis of 90°), and initial spin period of 1 ms, short 
enough to have no influence on the late-time timing proper- 
ties. The two models that would represent the two objects 
under discussion are presented in the following. 

IE 2259+586. 

This X-ray pulsar is characterized by a period P = 6.97 s and 
period derivative P = 4.8 x 10"'^ s/s (Gavriil & Kaspi 2002), 
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Fig. 4. — Evolution of timing, thermal and magnetic variables for some choices of initial magnetic field configuration that reproduce the observed properties 
of IE 2259+586 (solid line) and PSR J1814-1744 (dashed line). The spin-down ages of 230 kyr and 85 kyr are reached at real ages of 100 kyr and 60 kyr, 
respectively {top left panel). In the upper top right panel the evolution of the redshifted temperature at the pole (thick) and the equator (thin) is shown for both 
models. 



implying a timing age of t = 230 kyr. Its X-ray luminosity 
in the 2-10 keV band of 3.4 x lO"*"* erg/s (for a distance of 
4 kpc, Tian et al. 2010), much larger than its rotational en- 
ergy, prompted its classification as an AXP. The discovery of 
a major outburst in 2002 (Kaspi et al. 2003) further tight- 
ened the connection between this AXP and the SGRs, and 
established its 'magnetar-like' nature, despite the low inferred 
dipolar field of Bp = 5.9 x 10'^^ G. The quiescent (pre-outburst) 
spectrum of this object is characterized by the superposition 
of a thermal and a power-law component. The former, if fit- 
ted with a simple blackbody (BB), yields a BB temperature 
of about 0.4 keV and a BB radius of about 5.3 km. The pulse 
profile is double peaked (see Woods et al. 2004 for a summary 
of the pre-outburst spectral properties). The modulation level 
(or pulsed fraction, PF) was measured in two XMM obser- 
vations prior to the outburst (Woods et al. 2004); of interest 
here are the two softest bands (0.3-1.0 and 1.0-2.0 keV), dom- 
inated by the thermal component; the measured PF was found 
to be 0.169±0.015 in the 0.3-1.0 keV band and 0.195 ±0.006 
in the 1-2 keV band in the first of the two observations, and 
0.215 ±0.006 in the 0.3-1.0 keV band and 0.225 ±0.003 in 
the 1-2 keV band in the second observation (quoted errors are 
Icr; note that there is some variability between the measured 
values in the two observations). 

Fig. m (solid lines) shows an example of such evolution. 
The initial magnetic field has components Bp = 2.5 x 10'"* G 
and Btor = 2.5 x 10"' G. In the upper left panel we see how the 
timing age of the source increasingly departs from the real 
age. At the measured spin-down age of 230 kyr, the real age 
of the object is 100 kyr At this age the predicted luminosity 
of the object (in its thermal component, computed by integra- 
tion of the temperature profile over the whole surface of the 



star) is about 2 x 10 erg/s, its period about 6.7 s, its spin- 
down rate about 4.6 x lO"'-' s/s, and its Bp about 6.7 x 10'^^ G. 
Note that our realistic NS model described in Section 2 has 
a different radius and moment of inertia from the canonical 
values (10 km, lO**^ g cm^) usually assumed to infer Bp. Our 
goal is not to make a formal fit to the properties of the object, 
but simply to show typical examples of representative objects 
within a class. Nevertheless, the quantities in Fig. 1 are a very 
close match to the corresponding observed properties of this 
object. For the thermal quantities (luminosity, temperature), 
the comparison should however be made with care. The mea- 
sured luminosity is in fact subject to the uncertainty in dis- 
tance, and is also model-dependent. Our simulations also pre- 
dict the temperature profile during the pulsar evolution. We 
show the evolution of the surface temperature at the pole and 
the equator in the upper right panel of Fig. 4. A direct com- 
parison between our results and the observations can only be 
done at a qualitative level at this stage. In fact, while our pre- 
dictions specifically refer to the thermal luminosity, the mea- 
sured X-ray luminosity has both a thermal and a non-thermal 
(powerlaw, PL) component (e.g. Woods et al. 2004). The 
flux ratio PL/BB is a function of phase, varying between 0.3 
and 0.55 in the 0.5-7 keV energy band. If the non-thermal 
part is due to reprocessing of the thermal photons, then the 
comparison should be made with the total luminosity, but if 
the non-thermal photons are of different origin, then only the 
thermal component should be considered. 

In addition, the temperature inferred at infinity depends on 
both the compactness ratio M/R of the star, and the assumed 
model atmosphere, which greatly affects the spectrum. For 
the same NS compactness ratio (and hence same gravitational 
redshift), fits with Hydrogen atmospheres yield lower effec- 
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tive temperatures Tatm- The exact ratio rBB/7itm depends on 
a number of factors, such as atmosphere composition, mag- 
netic field strength, temperature, and, to a lesser extent, on 
column density and instrument type (see e.g. Suleimanov et 
al. 2011 for recent work using non-magnetic model atmo- 
spheres with varying composition). An extensive exploration 
of TfiB/Tatm for different values of the above parameters was 
made by D. Lloyd with detailed magnetized Hydrogen atmo- 
sphere models (Lloyd 2003). For a Z?-field strength of about 
10'^ G, a surface temperature of 1 - 2 x 10^ K, and a column 
density ~ 10^' cm"^, he found the ratio T^-R/T-^xm to be « 2 
(using the ACIS-1 spectra simulator). The observed spectrum 
of IE 2259H-586 has been fitted with a BBh-PL before the out- 
burst (Woods et al. 2004), and with both a BBh-PL and a 2BB 
after the outburst, as it approaches quiescence (Zhou et al. 
2008). In the observation closest to quiescence, Zhou et al. 
find that Tbb = 0.4 keV in the BBh-PL fit, and Tbb = 0.37 keV 
for the cooler component in the 2BB fit. Within these ob- 
servational uncertainties, and accounting for the temperature 
corrections due to a magnetized Hydrogen atmosphere, our 
theoretically predicted temperature would be roughly consis- 
tent with the BB temperature measured for this object in qui- 
escence. 

Further constraints on our predictions for the thermal vari- 
ables can be made through a comparison with the pulse profile 
and pulsed fraction of IE 2259+586. We used the predicted 
temperature distribution to compute the expected flux from 
the pulsar during its rotation, and hence calculate the pulse 
profile and the pulsed fractions in selected energy bands. For 
this calculation, we included the general relativistic effects of 
gravitational redshift and light deflection, following the for- 
malism of Page (1995), with the modifications by Perna & 
Gotthelf (2008). These allowed for the presence of a local 
anisotropic radiation pattern, f{5) oc cos 5" {5 angle between 
the normal to the surface and the outgoing photon direction, 
and n a 'beaming' parameter, the only free parameter that we 
used to find a match to the observed level of modulation), as 
suggested by realistic atmospheric models. Realistic models 
were limited to fields perpendicular to the star surface, (e.g. 
van Adelsberg & Lai, 2006), and hence not directly appli- 
cable to the spectrum from the entire NS surface. However, 
we still use them to gauge the extent of the radiation beam- 
ing. Using the models of Van Adelsberg & Lai (2006), we 
find that, for B 10'"^ G and T ~ 0.4 keV, the intensity for 
£■ ^ 1 keV drops as f{5) oc cos^ 5 up to 5 < 50 deg and is 
shallower at larger (5's. The effect of absorption to the source 
(parameterized by the column density Nh) was also included 
in our calculations, since it has been shown to be important 
for the computation of the PF in finite energy bands (Perna et 
al. 2000). 

We found the pulse profile to be double peaked, consistent 
with observations (Woods et al. 2004). We note that, for an or- 
thogonal rotator (and more generally for a large angle between 
the line of sight and the rotation axis), the "double peaked" 
nature of the pulse profile is a very robust prediction of our 
temperature profile, since it is symmetric with respect to the 
equator. For n = 1 .5, the pulsed fractions produced by the tem- 
perature profile of our simulations are 0.17 in the 0.3-1 keV 
band, and 0.22 in the 1-2 keV band. These values comfort- 
ably sit within the measured quiescent values for this source. 
However, we need to note that, given the BB plus powerlaw 
spectrum of this object, for the lowest energy band to be dom- 
inated by the thermal component, the powerlaw needs to have 
a cutoff in the soft X-ray band. This might be the case if hard 



photons are due to upscattering of the thermal ones. We did 
not attempt to model the pulse profile in a quantitative fash- 
ion, given the substantial contamination from the non-thermal 
component, for which no predictions of its local angular dis- 
tribution exist. A full spectral and timing analysis of the qui- 
escent spectrum (as in Perna & Gotthelf 2008 and Gotthelf et 
al. 2010) is beyond the scope of this paper, and could only be 
performed with a better handle of the physical origin of the 
powerlaw. 

PSRJ1814-1744 . 

This object has a measured period P = 3.975 s and period 
derivative P = 7.4 x 10"^^ s/s (Camilo et al. 2000), yielding 
a timing age of 85 kyr. The dashed lines in Fig. |4]show the 
evolution of a similar object. The initial magnetic field set up 
in this case is Bp = 1 .6 x 10'"* G and Bt™ = 8 x lO''* G. At the 
measured spin-down age of 85 kyr, the real age of the object 
is 60 kyr, with a predicted luminosity of about 4 x lO-'"' erg/s, 
5 times smaller than in the previous case. At the same age, 
its period is 4.4 s, its spin-down rate 8 x 10"'^ s/s, and its 
Bp about 9 X 10'-^ G. This is again in good qualitative agree- 
ment with the general properties of the object. Note that, for 
the X-ray luminosity, only upper limits exist. An analysis of 
ROSAT and ASCA archival data by Pivovaroff et al. (2000) 
yielded the values Lo.i-2.4keV < 6.3 x lO-'^tZ/lO kpc erg s"' 
from ROSAT L2-iokeV < 4.3 x 10^^c//10 kpc erg s"' from 
the ASCA band. Our predicted thermal (bolometric) luminos- 
ity is consistent with these limits. Given the lack of X-ray 
detection, no comparison can be made at this stage with the 
pulse profile. 

With these results, we have answered one of the puzzles: 
despite the second object under study has a similar dipolar 
field and is younger than the first example, its luminosity is 
substantially lower. Indeed, at other stages of evolution it can 
be even one order of magnitude smaller. The main reason 
for this difference, as we have discussed, is the larger initial 
toroidal field of the first example. 

Now we face the next problem: can the two objects have 
a markedly different behaviour in their outbursting proper- 
ties? To answer this question we have performed two sim- 
ulations as described in Sect. 3, following the evolution of 
the two objects for another 10"* yr but now using the code 
including the Hall term in the induction equation. For the 
same period of time, we found that the first object undergoes 
753 starquakes, while the second one only had 23 of these 
events. On average, one would expect that the first kind of 
source shows an outburst every 10-20 years, while the second 
object only displays a magnetar-like activity every 4-5 cen- 
turies. Our results are summarized in Fig.|5] In addition to the 
much lower event rate, outbursts of objects like PSR J1814- 
1744 are also less energetic than the average magnetar out- 
burst, thus making them harder to be observed. Therefore, 
the apparent dichotomy of the two objects with similar timing 
properties can easily be attributed to the slightly different ini- 
tial poloidal field and, more importantly, to the much stronger 
internal toroidal field of one of the NSs. 

5. SUMMARY 

In this work we have performed a comprehensive study of 
starquake s tatistics on the NS c rust, using the formalism de- 
veloped by Perna & PonsI (1201 1 ) for the coupled evolution of 
magnetic field and crust stresses during the lifetime of the NS. 
In particular, focusing on objects at age lO'*- 10^ yr, we have 
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Fig. 5. — Energy and waiting time distribution of the two cases representative of IE 2259+586(soIid) andPSR J1814-1744(dashes). The latter case is magnified 
a factor of 10 for comparison. The total number of outbursts corresponds to a period of ten thousand years. 



explored the dependence of the outburst frequency and en- 
ergetics on the initial strengths of the dipolar and toroidal 
fields. We have found that the outburst frequency is neg- 
ligible for initial poloidal fields below 10^"* G, even when 
the initial toroidal field is extremely strong, > 10^^ G . The 
rate of outbursts is an increasing function of Bp (at fixed 
Btor), and becomes significant (more than one per century) for 
Bp(t = 0) > 5 X 10^'' G. On the other hand, at fixed Bp, the ini- 
tial toroidal field has practically no influence on the starquake 
frequency for initial toroidal fields Btor ^5 ^p- This is because a 
toroidal field of the same order of magnitude of Bp is anyway 
rapidly formed on a HaU timescale. For initial Btor > Bp, the 
presence of the toroidal component affects the outburst rate, 
producing its rapid increase for initial Btor much stronger than 
Bp. 

The relative strengths of Bp and Btor at birth also influence 
the thermal quiescent luminosity of the NS in a similar fash- 
ion. For negligible Btor, the luminosity has a strong depen- 
dence on Bp(f = 0); it is about one order of magnitude larger 
for Bp ~ 10'^ G flian it is for Bp ~ 10'^ G. The influence of 
the toroidal field becomes important, as for the outburst fre- 
quency, only for Btor(? = 0) > Bp(t = 0). 

>From our simulations, we can conclude that there is a 
critical value of the initial dipolar field strength, Bp{t = 0) ^ 
3 - 5 X lO'"* G, above which some regular active periods may 
be expected. In this case, the intensity of the toroidal field 
at NS birth plays a crucial role in the ' astrophysical appear- 
ance' of those objects in their middle age. In fact, at around 
10^ yr, when Bp has roughly halved in strength to become 
~ 6-7 X 10^^ G, objects with neghgible Btoiit = 0) have an 
outburst rate below 1 per century, making them appear more 
like a quiet "high-B pulsar". On the other hand, with a larger 
toroidal field at birth, an object with similar Bp can have an 
outburst rate of a few years, making it appear as a 'magnetar' . 
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Correspondingly, due to the larger toroidal field, the thermal 
luminosity of the 'magnetar' is higher than that of the 'high-B 
pulsar' . Interestingly, observations have shown that middle- 
aged objects with measured Bp ^ 6-7 x lO'-' G can display 
radically different thermal and bursting properties. Specific 
examples are IE 2259+586 and PSR J1814-1744: while fliey 
both have Bp ~ 6 x 10'^ G, the former manifests itself as a 
magnetar, while the latter as a 'high-B pulsar'. Here we have 
shown how appropriate choices of Bp(f = 0) and Btor(f = 0) 
can naturally account for these differences, which no longer 
constitute a puzzle within the NS phenomenology. 

The apparent dichotomy between quiet, high B-field radio- 
pulsars and active magnetars is not real. There is a contin- 
uum of possibilities, and all types of sources can potentially 
show some unusual activity (e.g. the recent case of SGR 
0418+5729, Rea et al. 2010). We tend to classify them in 
classes according to a biased version of the anthropic princi- 
ple: if it has several active periods on a human lifetime, we 
name it a magnetar; if it does not show an outburst on the 
same timescale, we put it in the group of quiet, normal neu- 
tron stars. But this separation is not necessarily associated to 
fundamental differences in the nature of the sources or their 
internal physics. 
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